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Summary 

This paper examines mode generation and propagation characteristics of a 2-D cascade due to 
incident vortical disturbances using a time domain approach. Full nonlinear Euler equations are solved 
employing high order accurate spatial differencing and time marching techniques. The solutions show the 
generation and propagation of mode orders that are expected from theory. Single frequency excitations 
show linear response over a wide range of amplitudes. The response for multi-frequency excitations tend 
to become nonlinear due to interaction between frequencies and self interaction. 


Introduction 

The gust-cascade interaction problem has been studied extensively using semi-analytical and 
numerical approaches. These techniques mostly employ a frequency domain approach, examining one 
frequency at a time. With the availability of parallel processing algorithms, time domain approach has 
become feasible for gust-cascade interaction study. A time domain approach has the advantage that all 
harmonics of interest can be extracted from one solution and thus, may be able to mimic the real flow 
more closely. Also, linear/nonlinear regimes, self interaction, and multi-frequency interaction may be 
explored. In the nonlinear range, energy transfer between different frequencies (harmonics) occur, and 
such energy transfers are easily handled in a time domain approach. 

The 2-D cascade considered here is an unrolled section at a radial station of a modern high speed 
turbofan stator (fig. 1), for which flow and noise data are available. The incident gust is the periodic mean 
wake impinging on the stator vanes. The measured wake is represented using a Fourier series which 
includes only three harmonics of the blade passing frequency (BPF). The acoustic response of the cascade 
for this specified gust is studied employing a time domain approach. The full nonlinear Euler equations 
are solved employing high order accurate spatial differencing schemes, time marching schemes, and 
boundary conditions. 

In an earlier study (ref. 1), the acoustic response of the above 2-D cascade was examined for single 
and two frequency excitations. It was found that the acoustic response is linear for single frequency 
excitations. In the case of two frequency excitations (BPF+2BPF), the response was linear only when the 
BPF amplitude was small. For higher amplitudes of BPF, the response was nonlinear due to self 
interaction. The results also showed that the response was linear for equal excitation amplitudes of the 
two harmonics. When the amplitudes of the two harmonics are unequal and BPF amplitude is high, the 
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self interaction of BPF (which is cutoff for the geometry) influences the amplitude of the propagating 
2BPF mode in a significant way, both in the inflow and outflow regions. 

The 2-D gust-cascade interaction was one of the problems in the Fourth Computational Aeroacoustics 
Workshop on Benchmark Problems 2 . There, the acoustic response was sought for excitation at 
3 harmonics (BPF, 2BPF, and 3BPF). The BPF amplitude of excitation was specified such that acoustic 
response would be in the linear range. While BPF response should show an exponential decay (cutoff), 
the propagating mode amplitudes at 2BPF and 3BPF are of particular interest. At 2BPF, the 
circumferential mode order m = - 5 is a well cuton mode and propagates in the inflow and outflow regions 
(see fig. 2 for definition of these regions) without the influence of BPF. At 3BPF, the theory predicts 
propagation of two circumferential mode orders, m = 6 and m = -21. While mode m = 6 is well cuton, 
mode m = -21 is very close to cutoff. The flow solution obtained by the present time domain approach 
with the downstream boundary located at x = 1.5 (as specified in the problem, see fig. 2) showed that the 
m = - 21 mode propagates upstream as expected but at the outflow shows a slow decay and appear to 
propagate only near the end of the computational domain. This was examined further by obtaining a fine 
grid solution which also produced mode propagation characteristics as that of the coarse grid solution. 

Next, a grid stretched gradually in the x-direction (fig. 3) of the outflow region to x = 12.85 was used 
so that no reflected component, if any (from the downstream boundary) can influence the mode 
propagation characteristics, in particular that of the mode close to cutoff m = -21. The solution obtained 
with this stretched grid, showed that at 3BPF, m = - 21 mode also propagates in the inflow as well as 
outflow regions. 

In the present paper, the stretched grid is used to study the linear/nonlinear behavior of single and 
multi-frequency excitations. The acoustic responses obtained over a range of amplitudes of excitations are 
presented and discussed. 


2-D Cascade Problem 

A sketch of the gust-cascade problem is shown in figure 1 . The velocity triangle at the inflow shows 
the gust (representing the wake from the rotor blades) in the direction of relative flow velocity ( AT) at an 
angle (3 to the x-axis. The direction of tangential velocity (Qr) is in the positive y-direction as indicated 
resulting in the axial velocity (Too) making an angle a to the x-axis. The cascade has a gap-to-chord ratio 
of d/c = 2/3 with the inflow and outflow planes located at x = 3/2 c, where c is the vane chord, and d is the 
vane gap. The mean flow conditions at the inflow and outflow planes are given as: 

Inflow: Pi = 1, Ti = 1, and = 36° 

Outflow: p 0 l Pi = 0.92 

where P f and 7) are the normalized mean stagnation pressure and stagnation temperature respectively. a t is 
the mean flow angle and p 0 is the mean static pressure at the outflow plane. The inflow periodic wake 
disturbance is described at the inflow plane as: 

(jV,/) = |«1 cos(k y y- CO/) + a 2 cos^^y-©*)) + a 3 cos(3(£ y j>-©*))}^ 

p*Ov) =0 > Pg{y> t ) = ° 

% =cos(p)e I — sin(p)e J ,, p = 44° 


co 37i/4, k y = 1 Ik/ 9, a\ 0.005, a 2 = 0.003, a 3 = 0.0007 
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where co is the fundamental reduced frequency, k y is the transverse wavenumber, and a t s are the gust 
harmonic amplitudes. 

The frequency is normalized by the chord divided by the ambient speed of sound, wavenumber is 
normalized by the vane chord, and gust amplitudes are normalized by the ambient speed of sound. The 
harmonic amplitudes of excitation a x (at BPF), a 2 (at 2BPF) and a 3 (at 3BPF) have been varied to study 
the nonlinearities associated with the high amplitudes and multi-frequency excitations. 


Governing Equations 

The full nonlinear Euler equations governing the 2-D cascade flow are solved in the present 
investigation. In Cartesian coordinates these equations are written as: 

dQ dE dF n 
dt dx dy 


where 


P 


Q = 


p u 
pv 
E 


( 1 ) 


( 2 ) 


E = 


p u 

p u 2 + p 
p uv 

u(E + p ) 


pv 
p uv 

pv 2 + p 
v(E + p) 


(3) 


(4) 


and 


^ = (y-l)|£-lp( W 2 +v 2 )| (5) 

These equations are cast in generalized curvilinear coordinates employing the chain rule formulation as: 

3t dxd^ dxdv[ dyd^ dydr \ 
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Solution Procedure 


The full nonlinear time dependent Euler equations governing the 2-D cascade flow are solved 
employing a parallel computational aeroacoustic (CAA) code developed by Hixon et al. described in part 
in reference 3. Some details of the computational techniques are discussed in references 4 and 5. The flow 
equations written in chain-rule curvilinear form are solved using a prefactored sixth-order compact 
scheme for spatial differencing. The time marching uses a 2N Storage fourth-order nonlinear extension of 
Hu’s 5-6 Low Dissipation and Dispersion Runge-Kutta (LDDRK) scheme. 6,7 An explicit 10th order 
constant coefficient artificial dissipation 8 is used at every stage of the Runge-Kutta solver to dissipate 
unresolved waves. On the airfoils, the time derivative of the velocity normal to the wall is set to zero. 

At the inflow and outflow planes Giles nonreflecting boundary conditions 9 for the unsteady flow are 
implemented. The characteristics are defined as in reference 9 for the inflow and outflow boundaries. 
While fourth order boundary conditions are used at the inflow, second order boundary conditions are 
employed at the out flow (fourth order boundary condition at the outflow was found to produce numerical 
instabilities). While computing the derivatives of the characteristics, the gust components are subtracted 
at the inflow boundary. Once the time derivatives of the velocity components are computed, the gust 
components are added to them. Changes in the conserved variables as used in the code are then computed 
at the inflow and outflow boundaries. 

The specified mean flow boundary conditions have to be satisfied at the inflow and outflow 
boundaries. To apply the mean flow boundary condition (MFBC), the mean flow is directly computed at 
the boundaries as the flow evolves. The time derivative of the flow at the boundary is written as a sum of 
three components: 


(^)boundary ( (2 /)mFBC (Qt ) gust (^)non-reflective 

The method of evaluating (£?*)mfbc is discussed in detail in reference 10. This solution procedure has 
been successively used to study the gust response and the influence of steady loading. 11,12 


Computational Domain and Numerical Solution 

For the cascade configuration considered here the blade passing frequency is cutoff, at 2BPF the 
propagating mode is m= -5, and at 3BPF the propagating modes are m = 6 and m = - 21 (ref. 2). The grid 
used in this study is a stretched grid (fig. 3) to dissipate the outgoing waves at large distances in the 
downstream direction. The grid has 8 blocks with a total of 9506 grid points in each flow passage. Care is 
taken to increase the grid resolution in the leading edge, trailing edge, and wake regions. In the current 
solution procedure, 27 passages are solved and a periodic boundary condition in the y-direction is 
specified. With the specified gust amplitude, the solution is run until periodicity in pressure is achieved 
on the airfoil and on the inflow and outflow boundaries. The solution was run up to 235 periods so that 
change in circumferential mode amplitudes between successive periods is within one percent. Once the 
periodic solution is established, the solution is processed to obtain the acoustic response of the 
propagating modes at 2nd and 3rd harmonics. 


Joint Temporal-Spatial Transform 

The cascade unsteady flow solution as a function of time is known over a specified period. The 
acoustic pressure response is not only harmonic in time, but also periodic in the tangential direction, y. 
From the known acoustic pressure p as a function of (x, y, t), the acoustic mode amplitude P as a function 
of (x, m, n ), where m is the circumferential mode order and n is the harmonic order, is written 2 as 
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P(x,m,n ) 


2 




/=0 7=0 


where, Af is the time between samples, N Al is the number of samples over interval T 7 , frequency /= nIT ' 
and « = 0, 1, 2, . . ..(TV A/ - l)/2 is the harmonic number. TV is the number of points in the tangential direction, 
m is the spatial (circumferential) mode order, -TV/ 2 < m < (TV/2 - 1). 


Results and Discussions 

The cascade acoustic response for a vortical disturbance is studied using the time domain approach. 
The incident vortical disturbance includes three harmonics. The acoustic response of the propagating 
modes at 2nd and 3rd harmonics for single and multi-frequency excitations are examined. 


Propagating Modes 

For the cascade-gust interaction problem considered here, theory predicts propagation of m = -5 
mode at 2BPF and two modes in = 6 and in = -2 1 at 3BPF. However, in = -2 1 mode is very close to 
cutoff. The initial grid (fig. 2) used for the computation had the computational domain in the x-direction 
as specified in the benchmark problem, 2 from x = -1.5tox = 1.5. The solution obtained with this grid 
showed propagation of all the expected modes at the inflow. At the outflow, only m = - 5 and 6 modes 
propagated while m = - 21 mode was found to decay slowly. It was suspected that the small reflections 
from the downstream boundary with the current implementation of the nonreflecting boundary conditions 
may contribute to the observed behavior of the mode m = - 21 . To avoid those reflections, the grid was 
stretched in the downstream direction to x = 12.85 (fig. 3). The solution obtained with this grid for the 
excitation amplitudes specified by the benchmark problem showed that the mode m = - 21 also clearly 
propagated as expected from theory. The propagating mode amplitudes from the two solutions are 
compared in figure 4 for the inflow and outflow regions. At the inflow, the two solutions produce nearly 
similar amplitudes of the three propagating modes. At the outflow, m = - 5 and m = 6 mode amplitudes are 
the same for the two solutions. The effect of small reflections on the axial variation of 2BPF amplitude is 
not seen with the stretched grid solution and in = -2 1 mode clearly propagates in the downstream 
direction. 

Figure 5 shows the complex amplitude of the propagating modes in the inflow and outflow regions. A 
perfectly propagating wave would appear as a circle in such a plot. At 2BPF (fig. 5(a)), m = - 5 mode 
propagates well in the upstream and downstream directions. In figure 5(b), complex amplitudes of the two 
modes (m = 6 and -21) are shown for the inflow and outflow regions and they clearly show well 
propagating modes. Hence in the present study, the stretched grid is used to examine the linear/nonlinear 
behavior of single/multi frequency excitations. 

The pressure distributions on the airfoil surface, obtained with regular and stretched grids, for the 3rd 
harmonic excitation (a 3 = 0.0007) are shown in figure 5(c). The plots show substantial differences in 
pressure along the airfoil chord between the two solutions. While differences in magnitude are seen 
through the entire chord, the characteristics are different only near the leading edge region. Even here the 
leading edge peaks are nearly identical. But beyond the leading edge peak, loading difference occurs. It 
appears that this change in loading near the leading edge region is instrumental in changing the 
downstream propagation characteristics of m = -21 mode. 
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Single Frequency Excitation: Effect of Amplitude of Excitation 


Acoustic response characteristics of the cascade, for single frequency excitation is first explored. 
Figure 6(a) shows amplitudes of the propagating mode m = - 5 for 2BPF excitation in the upstream and 
downstream directions. When the excitation amplitude is increased 10 times to a 2 = 0.03, the 2BPF 
response also increases 10 times as shown in the figure. The dashed line represent the response for the 
excitation amplitude a 2 = 0.003 multiplied by a factor of 10. The linearity of the response for 2BPF 
excitation is clearly exhibited. 

For 3BPF excitation, a similar increase in excitation amplitude by a factor of 10 to 0.007 (from the 
original 0.0007) is shown to increase amplitudes of the propagating modes m = 6 (fig. 6(b)) and m = - 21 
by the same amount (fig. 6(c)) in inflow as well as outflow regions. Again, the dashed lines represent the 
response for the excitation amplitude of a 3 = 0.0007 multiplied by a factor of 10. The acoustic response of 
the cascade for 3BPF is clearly seen to be linear. 

Figure 7 shows a plot of excitation amplitude against the acoustic response for upstream and 
downstream propagating waves. One interesting feature to note is that in the downstream direction for 
equal amplitudes of excitation, 3BPF mode m = - 21 has higher response amplitudes than that of the 2BPF 
( m = -5) which was the dominant one for the original amplitudes of excitation specified in the problem. 


Multi-Frequency Excitation 

Next, the amplitudes of excitation in all three harmonics are specified to be 10 times the original 
amplitude. Figure 8 shows response amplitudes of the propagating modes for three-frequency excitation, 
ci\ = 0.05, a 2 = 0.03, and a 3 = 0.007. It is immediately seen that response for multi-frequency excitations 
are substantially different. At 2BPF, m = - 5 is well propagating at the inflow with a higher amplitude 
compared to 2BPF alone excitation. At the outflow, while the amplitude is much different, it appears to 
depart from linear behavior (fig. 8(a)). 

At 3BPF, in the flow region, m = 6 mode amplitude is smaller while m = - 21 mode amplitude is 
higher compared to 3BPF alone excitation. At the outflow an opposite trend is observed, i.e., the 
amplitude of m = 6 mode is higher and m = - 21 mode is lower than that of 3BPF only excitation, and 
appears highly nonlinear (fig. 8(b)). The modification of propagating mode amplitudes and the departure 
from linear behavior appear to stem from self interaction of BPF. 

To examine self interaction effects further, keeping the excitation amplitudes of 2BPF at 0.003 and 
3BPF at 0.0007, only the BPF amplitude of excitation is varied. The results of this study are shown in 
figure 9. It is seen that increase/decrease of the each propagating mode follows the pattern described 
above with reference to figure 8. But the departure from linearity starts to appear particularly at the 
outflow in all the propagating modes when the BPF amplitude of excitation is 0.02. The departure from 
well propagating modes is shown in figure 10, where complex amplitudes of the downstream propagating 
modes are shown. It is more clearly seen with respect to modes m = -5 and m = 6. 

When the amplitudes of excitation of the first and second harmonic are kept equal and increased 
gradually, the characteristics of the propagating modes are quite different as shown in figure 1 1. At 2BPF, 
the mode m = - 5 is well propagating and linear (fig. 1 1(a)) as would be expected since a x = a 2 . The 3BPF 
response amplitude undergoes complicated modifications due to self interaction and interaction between 
harmonics (frequencies). Also shown for comparison are the response amplitudes of single frequency 
excitation. Only upstream propagating m = - 21 mode remains linear and shows monotonic increase in 
response with increasing BPF+2BPF excitation amplitude. 

The linear behavior of 2BPF and departure from linear response of 3BPF are illustrated in figure 12, 
by the complex mode amplitudes of the propagating modes at the outflow region. At 2BPF mode m = - 5 
is clearly well propagating while at 3BPF, m = 6 shows first signs of departure from linear behavior of a 
well propagating mode. 
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Equal Amplitude Excitation 


When the amplitudes of excitation are equal, then the propagating modes are expected to remain 
linear and well propagating. This is illustrated for excitation amplitudes of a x = a 2 = a 3 in figure 13. Also 
shown are the response amplitudes for excitations with a 2 = a 3 and a x = 0 (dashed lines). The results 
clearly show that equal amplitude excitations result in nearly linear response and a x has little influence on 
the propagating modes. 


Nonlinear Effects: Loading and Wakes 

The factors that initiate the nonlinear responses are explored in this section. Figure 14(a) shows the 
typical loading on airfoil surface for three harmonics of excitations. Although the first harmonic loading 
is high, the BPF modes are cutoff for this configuration. The magnitude of loading decreases for higher 
harmonics. 

Next the 2BPF loading variations, for single frequency excitation, are examined. Figure 14(b) shows 
the loading for 2BPF excitation at a 2 = 0.003. Also shown are the 2BPF loading due to self interaction 
when excitation is at BPF and the amplitude is high. Loading for two amplitudes of excitation a x = 0.015 
and a x = 0.03 are shown in the figure. It is seen that the character (shape) of the loading due to self 
interaction is different from that of the 2BPF alone excitation. Also, the loading due to self interaction due 
to BPF excitation with a x = 0.03 is higher than that of 2BPF alone excitation. This explains at least in part 
the high amplitudes of acoustic response when the amplitude of BPF excitation is increased (fig. 9). 

Figure 14(c) shows the 2BPF loading for two frequency excitation (BPF+2BPF) compared with 
loading due to 2BPF only excitation. When the amplitudes of excitation are unequal, a x ^a 2 , the character 
of loading is different from the loading at amplitudes of excitation in the linear range. For equal 
amplitudes of excitation, a x = a 2 = 0.03, the character of the loading curve is similar to that of loading in 
the single frequency linear range excitation (at a 2 = 0.003). Hence, the change in character/shape of the 
loading curve is also an indication of departure from linear behavior. And, self interaction tends to change 
the nature of the loading curve. 

The downstream propagating waves tend to depart from linearity at amplitudes lower than upstream 
propagating ones. This is due to the wake development behind the vane trailing edge. The thickness of the 
wake increases with increasing amplitude of excitation as illustrated in figure 15, which shows the 
vorticity contours in the wake region of airfoil trailing edge for three amplitudes of excitation. 


Concluding Remarks 

A time domain approach was employed to study the gust cascade interaction problem. In this 
approach, full nonlinear Euler equations are solved using high order accurate spatial differencing and time 
marching techniques. Single frequency and multi-frequency excitations were considered to examine 
linear/nonlinear effects. The results of single frequency excitations show linear response over a wide 
range of amplitudes. Multi frequency excitations show that interaction between frequencies and self 
interaction modify the 2BPF and 3BPF responses. The higher harmonic responses become non monotonic 
and eventually exhibit nonlinear behavior. A complicated nonlinear response of 3BPF is observed when 
excitation amplitudes of BPF and 2BPF are equal. On the other hand, when the amplitudes of excitation at 
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the three harmonics are equal, nearly linear behavior is retained. The change in characteristics of airfoil 
loading and increasing thickness of the wake behind the airfoil trailing edge appear to contribute to the 
observed nonlinear behavior. Actual measurements may include such nonlinear effects and the present 
approach may help in understanding such features. However, linearized Euler analysis 13 has been shown 
to produce reasonably well the first two harmonics of rotor-stator interaction noise. The time domain 
approach will be used next to compute rotor-stator (3-D) interaction noise and compare with the 
experimental data. Such a study will clarify the influence of nonlinear effects on computed noise levels in 
annular cascades of interest. 
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Figure 4. Amplitude of propagating modes: regular and stretched grid solutions 
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Figure 5(a). Complex pressure amplitudes of propagating modes at 2BPF 




Figure 5(b). Complex pressure amplitudes of propagating modes at 3BPF 
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Figure 6(a). Single frequency excitation: linear behavior at 2BPF 



Figure 6(b). Single frequency excitation: linear behavior at 3BPF 
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Figure 7. Single frequency excitation: excitation amplitude and acoustic response of propagating modes 
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Figure 8(a). Multi-frequency excitation: Propagating mode amplitudes at 2BPF 






Figure 8(b). Multi-frequency excitation: Propagating modes at 3BPF 


NASA/CR — 2004-2 13110 


13 


Amplitude Amplitude _ p Amplitude 


0.001 





Axial distance Axial distance 




Figure 9. Multi-frequency excitation: propagating mode amplitudes 
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Figure 10. Multi-frequency excitation: Complex pressures of propagating modes, 

= 0.02, a 2 = 0.003, a 3 = 0.0007 
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Figure 11(a). Multi-frequency excitation: ai = a 2 = 0.01, 0.02, and « 3 = 0.0007: 2BPF mode amplitudes 
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Figure 11(b). Multi-frequency excitation: a x = a 2 = 0.01, 0.02, and «3 = 0.0007: 3BPF mode amplitudes 


NASA/CR — 2004-2 13110 


16 



Real P Real P 

Figure 12. Multi-frequency excitation: tti = « 2 = 0.01, a 3 = 0.0007: complex pressure amplitudes 
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Figure 13. Linear behavior of equal amplitude excitations 
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Figure 14(a). Typical harmonic loading on airfoil 



Figure 14(b). Single frequency excitation: 2nd harmonic loading 



Figure 14(c). Two frequency excitation: 2nd harmonic loading 
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Figure 15. Airfoil wake vorticity contours 
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